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Abstract 

Radiative transfer in a relativistic accretion disk wind is examined under the plane-parallel approxi- 
' mation in the fully special relativistic treatment. For an equilibrium flow, where the flow speed and the 

source function are constant, the emergent intensity is analytically obtained. In such an equilibrium flow 
■ the usual limb-darkening effect does not appear, since the source function is constant. Due to the Doppler 

and aberration effects associated with the relativistic motion of winds, however, the emergent intensity is 
strongly enhanced toward the flow direction. This is the relativistic peaking effect. We thus carefully treat 
and estimate the appearance of relativistic winds and jets, when we observe them in an arbitrary direction. 
Key words: accretion, accretion disks — astrophysical jets — gamma-ray bursts — radiative transfer 
' — relativity 
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^ ■ 1. Introduction 

H 

' Accretion disks are now widely believed to be energy sources in various active phenomena in the universe (see Kato 
Ci [ et al. 2007 for a review). Relating to energetic accretion disks, accretion disk winds have been extensively examined 
in connection with astrophysical jets and outflows: in bipolar outflows from young stellar objects (YSOs), in outflows 
from cataclysmic variables (CVs) and supersoft X-ray sources (SSXSs), in relativistic jets from microquasars (/iQSOs), 
^ I active galactic nuclei (AGNs), quasars (QSOs), and in gamma-ray bursts (GRBs). In particular, intense radiation 
— ' fields of luminous supercritical accretion disks may be responsible for relativistic jets from super-Eddington sources, 
00 . such as luminous juQSOs, GRS 1915+105 and SS 433, luminous QSOs, 3C 273, and energetic GRBs (see, e.g., Fukue 
CO ' 2004 for references). Furthermore, energetic emissions from relativistic jets have been examined, relating to, e.g., 
gamma-ray blazars and gamma-ray bursts (Dermer, Schlickciser 1993, 2002; Dermer 1998; Bottcher, Dermer 2002; 
Dermer et al. 2007). 

In such circumstances, radiative transfer in accretion disk winds as well as accretion disks becomes more and 
. more important. Radiative transfer in the standard disk has been investigated in relation to the structure of a 
static disk atmosphere and the spectral energy distribution from the disk surface (e.g., Meyer, Meyer-Hofmeister 1982; 
Cannizzo, Wheeler 1984). Furthermore, gray and non-gray models of accretion disks were constructed under numerical 
J. - treatments (Kfiz, Hubeny 1986; Shaviv, Wehrse 1986; Adam et al. 1988; Mincshige, Wood 1990; Ross et al. 1992; 
rS ; Shimura, Takahara 1993; Hubeny, Hubeny 1997, 1998; Hubeny et al. 2000, 2001; Davis et al. 2005; Hui et al. 2005) 
• and under analytical ones (Hubeny 1990; Artcmova ct al. 1996; Fukue, Akizuki 2006a). 

Radiative transfer in the accretion disk wind, on the other hand, has not been well considered both in the non- 
relativistic and relativistic regimes. For example, transformation properties of disk radiation fields in the proper 
frame of a relativistic jet were examined by, e.g., Dermer and Schlickciser (2002). In these earlier works, however, 
the radiation fields are set to be external sources, and the radiation transfer was not considered. Recently, radiative 
transfer in a moving disk atmosphere was firstly investigated in the subrelativistic regime (Fukue 2005a, 2006a), and 
in the relativistic regime (Fukue 2005b, 2006b; Fukue. Akizuki 2006b). In these studies, however, only the radiative 
moments were obtained under the moment formalism, and the specific intensity was not solved. Hence, in Fukue 
(2007) the specific intensity from an accretion disk wind was obtained in the subrelativistic regime, where the flow 
speed V is of the order of (v/c)^. 

In this paper, we thus extend the previous work. Namely, we examine radiative transfer in the relativistic accretion 
disk wind, which is assumed to blow off from the luminous disk in the vertical direction (plane-parallel approximation), 
with the relativistic speed up to the order of c. 

In the next section we describe the basic equations. In section 3, we show analytical solutions of the specific intensity. 
The final section is devoted to concluding remarks. 
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2. Relativistic Radiative Transfer Equation 



Let us suppose a luminous flat disk, deep inside which gravitational or nuclear energy is released via viscous heating 
or other processes. The radiation energy is transported in the vertical direction, and the disk gas. itself, also moves 
in the vertical direction as a disk wind due to the action of radiation pressure (i.e., plane-parallel approximation). For 
simplicity, in the present paper, the radiation field is considered to be sufficiently intense that both the gravitational 
field of, e.g., the central object and the gas pressure can be ignored. We also assume the gray approximation, where 
the opacities do not depend on the frequency. As for the order of the flow velocity v, we consider the fully special 
relativistic regime. 

The radiative transfer equations are given in several literatures (Chandrasekhar 1960; Mihalas 1970; Rybicki, 
Lightman 1979; Mihalas, Mihalas 1984; Shu 1991; Kato et al. 1998, 2007; Peraiah 2002; Castor 2004). The ra- 
diative transfer equation in the fully relativistic form is given in, e.g., the appendix E of Kato et al. (2007) in general 
and vertical forms. 

In a general form the radiative transfer equation in the inertial (fixed) frame is expressed as 
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Here, v is the flow velocity, c is the speed of light, and 7 (= 1/ -^/l — w^/c^) is the Lorentz factor. In the left-hand side 
the frequency-integrated specific intensity / and the direction cosine I are quantities measured in the inertial (fixed) 
frame. In the right-hand side, the mass density p, the frequency- integrated mass emissivity jo, the frequency- integrated 
mass absorption coefficient Kq^^, and the frequency-integrated mass scattering coefficient Kq^'^ are quantities measured 
in the comoving (fluid) frame, whereas the frequency-integrated radiation energy density E, the frequency-integrated 
radiative flux F, and the frequency-integrated radiation stress tensor P^^ are quantities measured in the inertial (flxed) 
frame. 

In the plane-parallel geometry with the vertical axis z and the direction cosine ^ (= cos0), the transfer equation is 
expressed as 
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where (3 (= v/c) is the normalized vertical speed, and F and P are the vertical component of the radiative flux and 
the radiation stress tensor measured in the inertial frame, respectively. 

For the convenience of readers, we shall show the full set of radiation hydrodynamical equations under the plane- 
parallel approximation, although we do not use all of them in this paper. 

For matter, the continuity equation, the equation of motion, and the energy equation become, respectively. 
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where u (= 7/3) is the vertical four velocity, J the mass-loss rate per unit area, ip the gravitational potential, e the 
internal energy per unit proper volume, p the gas pressure, and the internal heating. In the energy equation (5) 
the advection terms in the left-hand side are dropped under the present cold approximation. 

For radiation, the zeroth and flrst moment equations, and the closure relation become, respectively, 
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where f(T,P) is the variable Eddington factor, which is defined by / = Pco/Eco, Eco and Pco being the comoving 
quantities, and generally depends on the velocity and its gradient as well as the optical depth (Fukue 2006b, 2007b). 
Eliminating jo using the energy equation (5), the transfer equation (2) becomes 
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Introducing the optical depth defined by 
dT = ~ {k^^"" + kI''^) pdz, 
the transfer equation (9) finally becomes 
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Similarly, radiation hydrodynamical equations (4), (6), and (7), with the help of continuity equation (3) and the 
closure relation (8), become (cf. Fukue 2005b, 2006b) 
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Here, we dropped the gravitational and pressure forces. For such a radiation-dominated flow, where the gravitational 
and pressure forces are neglected, there are two integrals (Fukue 2005b, 2006b): 
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3. Analytical Solutions 



In order to solve the transfer equation (11) analytically, we suppose several assumptions; we assume that there is 
no internal heating (g"*" = 0) and the flow reaches the equilibrium state, where the flow speed is almost constant (/3 = 
const.). In this case, from equations (15) and (16), and the closure relation (8), the radiative flux F, the radiation 
stress tensor P, and the radiation energy density E are all constant. We further assume that the disk has a finite 
optical depth, and there exists a uniform isotropic source of intensity Iq at the disk equator of optical depth tq. 

Now, it is not so difficult to integrate the transfer equation (11). We first rewrite equation (11) symbolically as 
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where S'.^ + S'^ is the Doppler boosted source function and S'a + Ss is the non-boosted source function, both being 
independent of the optical depth under the present approximation. 

Under the above situations, we can formally integrate the transfer equation (11), similar to Fukue and Akizuki 
(2006) and Fukue (2007). After several partial integrations, we obtain both an outward intensity /(r,^,/3) (^ > 0) and 
an inward intensity /(r, — /i,/?) as 



/(r,M,/3) 



/(T,-/i,/?) 



7(1-/?/^) 



(r-ro) 



7(1 

+ /(ro,M)e ^^ 



(r-ro) 



7(1 



(20) 



(21) 



where I{To,fi) is the boundary value at the wind base on the luminous disk. These analytical solutions are reduced to 
those obtained in Fukue (2007) in the subrelativistic limit of 7 = 1. 

In general case with finite optical depth tq and uniform incident intensity Iq from the disk, the boundary value 
I{tq,^,P) of the outward intensity / consists of two parts: 

/(ro,/i,/3)=/o + /(To,-/i,/3), (22) 

where Iq is the uniform incident intensity and /(tq,— ^,/3) is the inward intensity from the backside of the disk beyond 
the midplane. Determining I{TQ,—fi,(3) from equation (21), we finally obtain the outward intensity as 
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Finally, the emergent intensity /(O,^,/?) emitted from the wind top becomes 
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In order to calculate the source functions, S'a and Ss, we consider two special cases below: a terminal case and an 
optically thin limit. 

3.1. Terminal Case 

When the flow speed is almost the equilibrium one and the radiation field is almost constant, then the values of the 
quantities of radiation fields are almost equal to those at the flow top. 

At the flow top of a moving photosphere at a relativistic speed, the usual boundary conditions for a static atmosphere 
is inadequate, as already pointed out in Fukue (2005b) Namely, the radiation fleld just above the wind top changes 
when the gas itself docs move upward, since the direction and intensity of radiation change due to the relativistic 
aberration and Doppler effect (cf. Kato et al. 1998, 2007; Fukue 2000). If a flat infinite plane with surface intensity Is 
in the comoving frame is not static, but moving upward at a speed Vs (= c/3s, and the corresponding Lorentz factor is 
7s), where the subscript s denotes the values at the surface, then, just above the surface, the radiation energy density 
Es, the radiative flux Fs, and the radiation pressure Pg measured in the incrtial frame become, respectively. 



cEs = 27r/s7, 
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Fig. 1. Normalized emergent intensity as a function of fi for several values of f3: (a) tq = 0.1 and (b) tq = 1 in the terminal case. 
The values of /3 are 0, 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, and 0.99; as the flow speed becomes large, the intensity at ^ = 1 is 
large. The dashed line is for the usual Milne-Eddinton solution for the plane-parallel case. 
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Fig. 2. Normalized emergent intensity as a function of fi for several values of tq in the case of /3 = 0.9 in the terminal case. The 
values of tq are 0.01, 0.1, and 1. The dashed line is for the usual Milne-Eddinton solution for the plane-parallel case. 
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In this case the non-boosted source functions arc calculated as 
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which depend on neither the flow speed nor direction cosine. Even if the non-boosted source function is constant, the 
emergent intensity (24) docs depend on the flow speed and the direction cosine. Examples of solutions arc shown in 
figures 1 and 2. 

In figure 1, the emergent intensity I{0,fi,P) normalized by Ig is shown for several values of /3 and tq as a function 
of 11. 

In the case of small tq (figure la), when the flow speed is small, the normalized emergent intensity is almost unity 
except for small fi direction, since the uniform source is seen except for small ^ direction, where the source function 
is seen. When the flow speed becomes large, however, the emergent intensity becomes remarkably anisotropic; it 
decreases in the edgeward direction, whereas it greatly increases in the poleward direction. This is the relativistic 
peaking effect, which originates from the relativistic Doppler effect and aberration. 

In the case of tq = 1 (figure lb), the emergent intensity for small (3 is reduced since the optical depth becomes large. 
However, the relativistic peaking effect becomes effective as the disk optical depth becomes large. 
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In figure 2, the normalized emergent intensity I{0,fj.,P) is shown for several values of tq in the case of /? : 
is seen in figure 2, the relativistic peaking effect strongly depends on the disk optical depth. 

3.2. Optically Thin Limit 

We next consider the flow at a constant speed using the closure relation (8) in the optically thin limit. 
Eliminating E by the closure relation (8), the source functions (18) and (19) are expressed as 
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is the scattering albedo. 

In the optically thin limit, we can easily obtain the quantities of radiation fields and the velocity- dependent Eddington 
factor (cf. Fukue 2006b, 2007b; Koizumi and Umemura 2007). That is to say, if there exists a uniform source of intensity 
Jo at the disk equator, the radiation energy density -E, the radiative flux F , and the radiation stress tensor P in the 
inertial (fixed) frame are respectively calculated as 



cE = 2tiIo, 

F = nIo, 
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Hence, from the closure relation (8) the velocity-dependent variable Eddington factor /(/3) is derived as 
' 3-3/3 + /32- 

Using these expressions, the source functions S'a and Ss are explicitly expressed as a function of P and /z: 
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The emergent intensity (24) does also depend both on the flow speed and the direction cosine. Examples of solutions 
are shown in figures 3 and 4. 

In figure 3, the emergent intensity /(0,/i,/3) normalized by Iq is shown for several values of f3 and tq as a function 
of fl. 

The qualitative behavior of the emergent intensity is similar to that of the terminal case. That is, in the case of 
small To (figure 3a), when the flow speed is small, the normalized emergent intensity is almost unity except for small 
fi direction, since the uniform source is seen except for small ^ direction, where the source function is seen. When the 
flow speed becomes large, however, the relativistic peaking effect becomes prominent. In the case of To = 1 (figure 3b), 
the emergent intensity for small (3 is reduced since the optical depth becomes large. However, the relativistic peaking 
effect becomes effective as the disk optical depth becomes large. 

In figure 4, the emergent intensity and the source functions normalized by Iq are shown for several values of (3. In 
figure 4, thin sohd curves represent Sg/A, S'a + Sg {A = 0.5), and <S'a/(l — A), from bottom to top, respectively, whereas 
thick solid curves denote I{0,n,/3) {A = 0.5). 

As is seen in figure 4 and equation (36), the source function S'a does depend only on the flow speed. This is just 
because, in the present dcflnition, the source function Sa is proportional to the radiation energy density E^o in the 
comoving frame, since cEco = cE- 2F/3 + f3^cP. On the other hand, the source function Ss depends on the direction 
cosine as well as the flow speed. It is small both in the forward and backward directions, while it becomes the maximum 
in the direction at = /3. 
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Fig. 3. Normalized emergent intensity as a function of for several values of /3: (a) ro = 0.1 and (b) ro = 1 in the optically thin 
limit. The values of /3 are 0, 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, and 0.99: as the flow speed becomes large, the intensity at ^ = 1 
is large. The dashed line is for the usual Milne-Eddinton solution for the plane-parallel case. 
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Fig. 4. Normalized emergent intensity and source functions as a function of fi for several values of /3 in the case of ro = 0.1 in the 
optically thin limit. The values of /3 are (a) 0.5 and (b) 0.9. Thin solid curves represent Sb/A, Ss, + Sb {A = 0.5), and Sa/ (1 — A), from 
bottom to top, respectively, whereas thick solid curves denote 7(0, ^,/3) {A = 0.5). The dashed line is for the usual Milne-Eddinton 
solution for the plane-parallel case. 



4. Concluding Remarks 

In this paper we have examined the radiative transfer problem in an accretion disk wind under the plane-parallel 
approximation in the fully special relativistic regime. For an equilibrium flow, where the radiative quantities and 
source function arc constant, we analytically obtain the specific and emergent intensities. We found that the emergent 
intensity depends on the flow speed as well as the direction cosine, and exhibits a relativistic peaking effect. As a 
result, a wind luminosity would be overestimated by a polc-on observer and underestimated by an edge-on observer, 
when we observe an accretion disk wind (cf. Sumitomo et al. 2007; Nishiyama ct al. 2007). 

It should be noted that the apparent optical depth in the relativistically moving media. Abramowicz et al. (1991) 
pointed out that the optical depth in the relativistic flow decreases as 7(1 — (3^)t toward the downstream direction, 
due to the Doppler and aberration effects. Inspecting equation (11) or solution (23), we find that in the present case 
the optical depth r is apparently replaced by 7(1 — (3^)t. This is just consistent with the results by Abramowicz et 
aL (1991). 

In this paper we only examined the frequency-integrated intensity under the plane-parallel and gray approximations. 
It should be briefly remarked on the frequency dependence; i.e., the frequency-dependent intensity 1^. As long as the 
opacity is gray, the spectral transformation is determined by the relativistic effect in the present situation. Namely, 
for the frequency-integrated intensity the relativistic invariant is I /v'^^ and this effect appears in solution (23) and 
other equations as a factor of [7(1 — /?^)]^. Since the relativistic invariant for the frequency-dependent intensity is 
ly/v^, the corresponding factor should be changed as [7(1 — /3^)]'^. On the other hand, the relativistic modiflcation 
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in the optical depth discussed above is not changed in the frequency-dependent case. Hence, the present results may 
be valid for the frequency-dependent emergent intensity; e.g., the incident intensity I^q would be boosted in the polar 
direction, according to solution (23), but [7(1 — is replaced by [7(1 — /3^)]'^. 

The present study can be applied to energetic jets in, e.g., gamma-ray blazars and gamma-ray bursts. The effect of 
the relativistic jets on the emergent spectrum has been studied in several literatures (Dermer, Schlickeiser 1993, 2002; 
Dermer 1998; Bottcher, Dermer 2002; Dermer et al. 2007), In these earlier works, however, the radiation transfer was 
not treated at all. Hence, the present approach may be very usefuU in these fields of active phenomena. 

The radiative transfer problem investigated in the present paper must be quite fundamental problems for accretion 
disk physics and astrophysical jet formation. In onder to demonstrate the existence of the relativistic peaking effect, 
we have imposed various assumptions, including a constant flow speed, gray approximation, no heating source, and 
so on. By relaxing these assumptions and integrating the relativistic transfer equation numerically, we could obtain 
the emergent intensity and spectra more quantitatively. These are left as future works. 

The author would like to thank S. Kato, S. Mineshige, M. Umemura, T. Koizumi, and C. Akizuki for enlightening 
and stimulating discussions. This work has been supported in part by a Grant-in- Aid for Scientific Research (18540240 
J.F.) of the Ministry of Education, Culture, Sports, Science and Technology. 
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